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Abstract. We aim to generate virtual commuting networks in the French 
rural regions in order to study the dynamics of their municipalities. Since 
we have to model small commuting flows between municipalities with a 
few hundreds or thousands inhabitants, we opt for a stochastic model 
presented by Gargiulo et al. (2012). It reproduces the various possible 
complete networks using an iterative process, stochastically choosing a 
workplace in the region for each commuter living in the municipality of 
a region. The choice is made considering the job offers in each munic- 
ipality of the region and the distance to all the possible destinations. 
This paper presents how to adapt and implement this model to generate 
French regions commuting networks between municipalities. We address 
three different questions: How to generate a reliable virtual commuting 
network for a region highly dependant of other regions for the satisfac- 
tion of its resident's demand for employment? What about a convenient 
deterrence function? How to calibrate the model when detailed data is 
not available? We answer proposing an extended job search geographical 
base for commuters living in the municipalities, we compare two different 
deterrence functions and we show that the parameter is a constant for 
network linking French municipalities. 



1 Introduction 

Some rural areas have an increasing population while others continue to suffer 
from depopulation in Europe (Johansson and Rauhut, 2007; Champcticr, 2000). 
This is what the European project PRIMA^ has been trying to understand 
studying the dynamics of a commuting network of virtual rural municipaties 
through microsimulation (Huct ct al., 2012). In this framework, as in many 
studies through microsimulations or agent-based simulations, we need genera- 
tion models capable of building reliable virtual commuting networks. Indeed, 
new economic theories assume local positive dynamics can be explained by im- 
plicit geographical money transfers made by commuters or retired people (see 

^ PRototypical policy Impacts on Multifunctional Activities in rural municipalities - 
EU 7th Framework Research Programme; 2008-2011; https:/ /prima. cemagref.fr/the- 
project 



for example Davczies (2009)). It is thus necessary to have virtual commuting 
networks of individuals for the micro simulation approaches increasingly used 
(Birkin and Wu, 2012; Parker et al., 2003; Bousquct and Page, 2004) to study 
the economic and land use dynamics. 

We are interested in understanding the dynamics of the French rural munic- 
ipalities. 95% of them have at most 3000 inhabitants. This means most of the 
commuting flows we want to study are weak, with a spatial distribution highly 
decided by chance. It's why we opt for a stochastic model proposed recently by 
(Gargiulo et al., 2012). Moreover, we want to consider the commuting network 
at different dates. Detailed data regarding flows between couples of municipal- 
ities are available in France only for 1999. At the other dates, the only reliable 
data is aggregated data, which for each municipality describes how many people 
go to work outside of the municipality and how many come from the outside 
of the municipality to work in without precisions about the various places of 
work or the various municipalities of residence. Then we also choose the model 
of Gargiulo et al. (2012) for its ability to generate a population of individuals 
on a commuting network, starting from this data. This model reproduces the 
complete network using an iterative process stochastically choosing a workplace 
in the region for each commuter living in the municipality of the region. The 
choice is made considering the job offers in each municipality of the region and 
the distance to all the possible destinations. It differs from the classical gener- 
ation models presented in Ortuzar and Willumsen (2011) since it is a discrete 
choice model where the individual decision function is inspired from the gravity 
law model which is not used usually at an individual level (Haynes and Fother- 
ingham, 1984; Ortuzar and Willumsen, 2011; Barthclcmy, 2011). Moreover, the 
model ensures that for every municipality the virtual total numbers of com- 
muters coming in and going out are the same as the ones of the data. This paper 
presents how to adapt and implement this model to generate French regions 
commuting networks between municipalities. This implementation has forced us 
to address three different questions: How to generate a reliable virtual commut- 
ing network for a region highly dependant of other regions to satisfy the need for 
jobs of people living in their municipalities? What about a convenient deterrence 
function? How to calibrate the model when detailed data is not available. 

A first problem we have to solve is that our French regions are not islands as 
for example in De Montis et al. (2007, 2010). Indeed, some of the inhabitants, 
especially those living close to the border of the region, are likely to work in 
municipalities located outside the region of residence. This part, especially if it 
is significant, make the generated network largely false if we only consider that 
people living in the region also work in the region. A way to solve this problem 
is to generate the commuting network only for people living and working in 
the region. However, this means the modeller has to know the quantity and 
the place of residence of individuals who work outside the region and live in 
the region. The data giving this detail is very rare, and so is being able to 
benefit from an expertise that could give this knowledge. Then, we address this 
issue by extending the job search geographical base for commuters living in the 



municipalities to a sufficiently large number of municipalities located outside the 
region of residence. We compare the model without municipalities from outside 
(called only outside later in this paper) and the model with outside on 23 French 
regions and conclude about the quality of our solution. 

The second problem relates to the form of the deterrence function which rules 
the impact of the distance on the choice of the place of work relatively to the 
quantity of job offers. The initial work done by Gargiulo ct al. (2012) proposes 
to use a power law. However, Barthclemy (2011) says the form of the deterrence 
function varies a lot, sometimes it can be inspired from an exponential function 
as in Balcan ct al. (2009) or from a power law one as in Viboud ct al. (2006). 
To choose the much more convenient deterrence function, we have compared the 
quality of generated networks for 34 French regions obtained on the one hand 
with the exponential law, and with the power law on the other hand. We showed 
that we obtained better results with the exponential law. 

The last, but not the least problem to solve is the one related to calibration. 
The generation model, as most of the currently used commuting network gener- 
ation models, has one parameter to calibrate. This parameter rules the impact 
of the distance in the individual decision regarding its place of work relatively to 
the quantity of job offers. The only available distance we can use is the Euclidian 
distance. We have detailed data on the commuting network for the year 1999 
which can be used for calibration, but it is not the case for early years or more 
recent ones. It may be possible to assume the parameter value does not change 
over time but we know the transportation network can evolve greatly at the local 
level to reduce the time distance while we can't capture such a change with the 
Euclidian distance. The solution was finally easy to manage. Using 34 French 
regions we show that every French region can be generating using a constant 
value for the parameter. Then, we assume that the parameter value is constant 
over time and space. 

2 Material and methods 

2.1 The French regions and the data from the French statistical 
office 

A complete description of the regions from which the network has been gener- 
ated is provided in the Table 4. These regions have been randomly chosen for 
their diversity in terms of number of municipalities and commuters, and surface 
areas. Some correspond to a region while others are closer to the county (called 
" departement" in French). 

The French Statistical Office (called IN SEE) collects information about 
where individuals live and where they work. From this collected data, the Mau- 
rice Halbwachs Center or the IN SEE make available for every researcher the 
following data: 

— in 1999, data about the numbers of individuals commuting from location i 
to location j for every municipality of a region; 



— in 1990 and 2006, the total number of commuters, the total job offers and the 
total number of workers living in for every municipality; these data allows 
to compute for each municipality the number of commuters coming to work 
in. 

The Lambert coordinates for each municipality are easy to find on inter- 
net. They allow to compute the Euclidian distance between every municipality 
couple. 

We start from these data sets for our implementation work of the model 
presented in the next section. 

2.2 The model of Gargiulo et al. (2012) 

Consider a region composed by n municipalities. We can model the observed 
commuting network starting from the matrix: R € M„xn(N) where Rij is the 
number of commuters from the municipality i (in the region) to the municipality 
j (in the region). This matrix represents the ligth grey origin-destination table 
presented in Table 1. 

The inputs of the algorithm are: 

— D = {dij)i<ij<n the Euclidean distance matrix between municipalities 

— Ij the number of in-commuters from the region to the municipality j of the 
region, 1 < ^ < n (i.e. the number of individuals living in the region in a 
municipality i {i ^ j) and working in the municipality j). 

— Oi the number of out-commuters from the municipality i of the region to 
the region, 1 < i < n (i.e. the number of individuals working in the region 
in a municipality j [j ^ i) and living in the municipality i). 

Ik and Ok can be respectively assimilated to the job offers for the employed of 
the region and the job demand of the employed of the region for the municipality 
k, 1 < k < n. The algorithm starts with: 

n 

and 

n 

0,:=^i?., (2) 

The purpose of the model is to generate the ligth grey origin-destination 
subtable of the region described in Table 1. To do this it generates the matrix 
S € M„xn(N) where Sij is the number of commuters from the municipality i 
(in the region) to the municipality j (in the region). It's important to note that 
Sij = a i — j. The algorithm assigns at each individual a place of work with 
a probability based on the distance from its place of residence to every possible 
places of work and their corresponding job offer. The number of in-commuters 



of the municipality j and the number of out-commuters of the municipahty i 
decreases each time an individuals living in i is assigned the municipality j as 
workplace. We stop the algorithm when all the out-commuters have a place of 
work. The algorithm is described in Algorithm 2.1 with m = n. 



Table 1. Origin-destination table for the region; The light grey table represents the 
commuters living (place of residence RP) and working (place of work WP) in the region 
for each municipality of the region; The dark grey line represents the number of out- 
commuters from a municipality of the region to the region for each municipality of the 
region (i.e. the row totals of the light grey table); The dark grey column represents 
the number of in-commuters from the region to a municipality of the region for each 
municipality of the region (i.e. the column totals of the light grey table). 
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Algorithm 2.1 commuting generation model 

Input : D G M„xm(K), / £ N™, O £ N", /3 G R+ 
Output : 5' G M„xm(N) 
S^j ^ 

while 0» > do 

Simulate i ~ where yl = {fc|fc G | [1, n]|, Ofc / 0} 
Simulate j from |[l,m.]| with a probability: 



Sij i Sij -\- 1 

^ - 1 
end while 
return S 



Gargiulo et al. (2012) uses a deterrence function f{dij,f3) with a power law 
shape: 

/(d„-,/3) = d-.^ l<i,j<n . (3) 



3 Statistical tools 



This section presents the tools used to caUbrate the model and to compare 
various implementation choices. 



3.1 Calibration of the (3 value. 

We used the same method as in Gargiulo ct al. (2012) to calibrate the f3 value. 
We calibrate j3 in order to minimize the average Kolmogorov-Smirnov distance 
between the simulated commuting distance distribution and the one building 
from the observed data. For the basic model we compute the commuting distance 
distribution with the commuting distance of the individuals who are commuting 
from the region to the region. For the model with the outside we compute the 
commuting distance distribution with the commuting distance of the individuals 
who are commuting from the region to the region and also to the outside. 

As the model of Gargiulo et al. (2012) is stochastic, the final calibration value 
we consider is the average /3 value over 10 replications of the generation process. 



3.2 An indicator to assess the change. 

We need an indicator to compare the simulated commuting network and the 
observed commuting network. Let R e M„jxn2(^) ^ commuting network when 
Rij is the number of commuters from the municipality i to the municipality j. 
Let S e M„jxn2(N) another commuting network of the same municipalities. We 
can compute the number of common commuters between R and S (Eq. 4) and 
the number of commuters in R (Eq. 5): 

ni 112 

A^cc„, x„, {s, R) = J2J1 ' -^^j) (^) 
i=l J=l 

i=l 3 = 1 

From (Eq. 4) and (Eq. 5) we compute the S0rensen similarity index (S0rcnsen, 
1948). This index makes sense since it corresponds to the common part of com- 
muters between R and S. Then we call it common part of commuters {CPC) 
(Eq. 6): 

C /^G„, {b,R)- ^^^^^^^ ^ ^^^^ (6) 

It has been chosen for its intuitive explanatory power: it's a similarity coef- 
ficient which gives the likeness degree between two networks. It is ranging from 
a value of 0, for which there aren't any commuters flows in common in the two 
networks, to a value of 1 when all commuters flows are exactly identical in the 
two networks. 



4 Generating French regions commuting network at 
municipality level 

4.1 How to cope with regions which are not island or with the lack 
of detailed data? 

A commuting network is defined by an origin-destination table (light grey table 
in Table 2). At the regional level, it means that we need to know, for each 
municipality of residence of the region and for each municipality of employment 
of the region, the value of the flow of commuters going from one to another. 
This kind of data is not always provided by the statistical offices and usually 
the datasets are aggregated: only the total number of out-commuters and in- 
commuters for each municipality is available for each municipality (dark grey 
row and colum in Table 2). To apply the model and define the commuting 
network, unless we are on a really isolated region^, we should need to find a way 
to isolate from the total number of in(out)-commuters (dark grey row and colum 
in Table 2) the fraction strictly relating to the region (light grey table in Table 
2). This is actually not a simple task. 

Moreover, even if we are able to isolate these parts, it remains a problem 
due to the border effect. Indeed, if we consider only the region, we risk to make 
an error in the reconstruction of the network of the municipalities close to the 
border of the region. The higher the proportion of individuals working outside 
of the region is, the higher the error will be. 

To go further, we propose to change the inputs of the algorithm. Instead 
of only considering the regional municipalities as possible places of work, we 
also consider an outside of the region. The outside represents the surroundings 
of the studied area. The following part describes how to consider this outside 
practically. 

A new extended to outside job search base. We implement the model, 
taking or not into account an outside, to generate 23 various French regions 
(see the 23rd regions in the table 4). Their outside is composed of the set of 
municipalities of their neighbouring '"departements"'. 

We consider the outside of the region composed by to — rt municipalities, 
where n is the number of municipalities of the region. The inputs arc the directly 
available aggregated data at the municipal level: 

— D — (dij) i<i<n the Euclidean distance matrix between the municipalities 

l<j<m 

both in the same region and in the outside 

— {Ij)i<j<rn the total number of in-commuters of the municipality j of the 
region and outside of it (i.e. the number of individuals working in the mu- 
nicipality j of the region or the outside and living in another municipality). 

~ (Oi)i<i<„ the total number of out-commuters of the municipality i of the 
region only (i.e. the number of individuals living in the municipality i of the 
region and working inan other municipality). 

^ an island for example, in this case grey rows and colums in Table 2 would not exist 



The purpose of the algorithm with introduction of the outside is to generate 
the origin-destination table (ligth grey and grey subtable in Table 2). To do 
this we use the algorithm presented in Algorithm 2.1 to simulate the Table 3. 
Now, we can easily obtain by difference the Table 2 with the total number of 
in-commuters {Ij)i<j<n, the total number of out-commuters (Oi)i<i<„ and the 
light grey table of the Table 3. 

We obtain a matricial representation of the origin-destination table pre- 
sented in the ligth grey and grey subtable in Table 2, the simulated matrix 
5 € M 

(n+i)x(ri+i)(I^) where Sij is: 

— the number of commuters from the municipality i (in the region) to the 
municipality j (in the region) ii i,j ^ n + 1; 

— the number of commuters from the outside to the municipality j (in the 
region) if i = n + 1 and j n + 1; 

— the number of commuters from the municipality i to the outside if i ^ n + 1 
and j = n + 1. 



Table 2. Origin-destination table; The light grey table represents the commuters living 
and working in the region for each municipality of the region; The grey column represent 
the out-commuters living in the region and working outside (Out.) for each municipality 
of the region; The grey line represents the in-commuters working in the region and 
living outside (Out.) for each municipality of the region; The dark grey line(column) 
represents the total number of out(in)-commuters for each municipality of the region. 
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Comparison of the two models: Assessing the impact of the outside. 

We assess the impact of the outside by a comparison between the generations of 
the network for 23 French regions with and without the outside. The generation 
is made at the municipality scale with a power law detterence function. 

The inputs of the case without-outside are built from detailed data while the 
inputs for the with-outside case are directly the aggregated data (the total mu- 
nicipal number of in and out-commuters). Both implementations are compared 
through their CPC values for each region. We replicate ten times the generation 



Table 3. Origin-destination table from the region to the region and the outside; The 
light grey table represents the commuters living (place of residence RP) and working 
(place of work WP) in the region for each municipality of the region; The grey table 
represents the commuters living (place of residence RP) in the region and working 
(place of work WP) outside of the region. 
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for each region and compute our indicator on each replicate. In all the pre- 
sented figures, the indicator is averaging on 10 replications. The variation of the 
indicator over the replications is very low, 1.02% of the average at most. Conse- 
quently, it is not represented on the figures. Fig. 1 presents the common part of 
commuters CPCnxniS, R) between the simulated network S and the observed 
network R obtained with the regional job search base (square) and obtained with 
a job search base comprising the region and its outside (triangle). It's important 
to note that for the implementation without outside 5* G M„xn(N) while for the 
implementation with outside 5* € M(„_|_i)x(n+i)(N). In order to compare the two 
models we just consider the regional network (commuters from the region to the 
region). Indeed, in the without outside case NCnx7i{S) = NCnxn{R) but it's 
not necessarily true for the with outside case. 

Fig. 1 shows that the two job search bases give results which are not really 
different. Thus, introducing the outside solves the problem linked to the lack of 
detailed data without changing the quality of the resulted simulated network. 
Indeed we have to keep in mind that the inputs of the with-outside case does 
not require to have detailed data on the contrary to the without outside case. 



□ 

|o.s- 




□ 



Fig. 1. Average CPC for 23 regions. 
The squares represent the basic model; 
The triangles represent the model with 
outside. 



4.2 Choosing a shape for the deterrence function 



The second problem relates to the form of the deterrence function which rules 
the impact of the distance on the choice of the place of work relatively to the 
quantity of job offers. The initial work done by Gargiulo et al. (2012) proposes 
to use a power law. However, Barthclcmy (2011) says the form of the deterrence 
function varies a lot, sometimes it can be inspired from an exponential function 
as in Balcan ct al. (2009) or from a power law one as in Viboud ct al. (2006). To 
choose the much more convenient deterrence function, we compare the quality of 
generated networks for 34 French regions obtained with the model with outside 
on the one hand with the exponential law, and with the power law on the other 
hand. 

A deterrence function following an exponential law is introduced: 

f{d,j,l3) = e"'^*^ 1 < i < n and 1 < j < m . (7) 

To compare the two deterrence functions, we have generated the networks of 
34 various French regions (see table 4 for details) replicating ten times for each 
region. The networks were generated with a job search base for the algorihtm 
considering the outside. 

As an example, the Fig. 2 shows we obtained a better estimation of the 
Auvergne commuting distance distribution with the exponential law. 



Fig. 2. Density of the Auvergne com- 
muting distance distribution; the solid 
line represents the observed commuting 
distance distribution; the dotted line 
represents the commuting distance dis- 
tribution obtained with the calibrated 
model with a job search base com- 
prising the outside and the exponen- 
tial law; the dashed line represents the 
commuting distance distribution ob- 
tained with a job search base compris- 
ing the outside and the power law. The 
two simulated commuting distance dis- 
tribution are computed for one replica- 
tion each. 



More systematically, we plot for the two different deterrence functions, ex- 
ponential law and power law, the average on the replications of the common 
part of commuters CPC(^n+i)y.{n+i){S,R) in the Fig. 3. It clearly shows that the 
average proportion of common commuters is always better with the exponential 
law represented by the squares. 
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Fig. 3. Average CPC for the power shape (triangle) and the exponential shape (square) 
for 34 french regions. 



4.3 Calibrating the model for French regions 

The last difficulty to solve is about the calibration process which requires until 
now to have detailed data to be accurate. 

Fig. 4 shows the calibrated /3 values for each of the 34 French regions. We can 
see these values weakly vary from about 1.7 • 10^'' to 2.4 • 10^* with an average 
(3 valued (C = 1.94 • 10""*) corresponding to the dark line. 

Then we hypothesize that it is possible to directly calibrate the algorithm to 
generate the 34 French regions, simply using a constant equal to C. To study 
the influence of this approximation on the common part of commuters we have 
computed the CPC with C as the parameter value for the 34 regions. We observe 
in the Fig. 5 that the influence of the /3's approximation on the CPC is very 
weak. We note that the average CPC obtained with C is, for some regions, 
higher than the CPC obtained with the not averaging /3 value. It's possible that 
the common part of commuters are better with another /3 value because it's not 
the calibration criterion. 

We don't need to study the influence of the /3's approximation on the cal- 
ibration criterion. Indeed from the studies made in Gargiulo et al. (2012), we 
know the CPC and the calibration criterion are highly correlated. The CPC 
and the calibration criterion have the same evolution in terms of /?. The /? value 
for the minimization of the Kolmogorov-Smirnov distance is very close to the 
one obtained for the maximization of the CPC (see the figure 7 in Gargiulo 
et al. (2012) which perfectly illustrates this relation). Then, as CPC values re- 
main quasi identical with (3=C or with /3 valued from the calibration process 



presented in 3.1, the quality of the approximation of the caHbration criterion, 
i.e. the commuting distance distribution, remains the same. 
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Fig. 4. The circle represents the 
average calibrated /3 values for ten 
replications (The confident inter- 
val is composed of the minimum 
and the maximum) for each re- 
gions; the line represents the av- 
erage P value for the 34 regions. 
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Fig. 5. Common part of commuters for the 34 regions; The squares represents the 
average CPC (10 replications) obtained with the calibrated j3 value; The triangles 
represents the average CPC (10 replications) obtained with the estimated /3 values 
(average /3 value over the 34 calibrated j3 values). 



5 Discussion and conclusion 

To study the rural areas dynamics by microsimulation, we need virtual commut- 
ing networks linking individuals living in the municipalities of various French 
regions. As the studied scale is very low, we have small flows and decided to opt 



for a stochastic generation algorithm. The one recently proposed by (Gargiulo 
et al., 2012) is relevant for our problem. Starting from this model, we implement 
the commuting network of 34 different French regions. The implementation work 
leads us to solve three practical problems. 

The first problem we have to solve is that our French regions are not islands. 
Indeed, some of the inhabitants, especially those living close to the border of the 
region, are likely to work in municipalities located outside the region of residence. 
However the classical approaches to generate commuting network consider only 
residents of the region working in the region. That is also the case of ours. The 
data giving details, or the knowledge, allowing the modeller to suppress people 
living in the region and working outside is hard to obtain. Then, we address 
this issue by extending the job search geographical base for commuters living 
in the municipalities to a sufficiently large number of municipalities located 
outside the region of residence. We compare the model without municipalities 
from outside and the model with outside on 23 French regions. We conclude 
about the relevance of our solution which keeps the value of our quality indicator 
identical while it does not obliged to know about people who don't work in the 
region and permit to generate networks only starting from aggregated data. 

(Gargiulo et al., 2012) model is based on the gravity law. Then, the sec- 
ond problem relates to the deterrence function which is more a power law or an 
exponential law depending on the study. Moreover, as the empirical studies com- 
paring the generated networks to "real" data are very rare (Barthclcmy, 2011), 
no one knows about the better shape. To choose the much more convenient one 
for our French regions, we have compared the quality of generated networks for 
34 regions obtained on the one hand with the exponential law, and with the 
power law on the other hand. We showed that we obtained better results with 
the exponential law whatever the region is. Indeed our 34 regions vary a lot in 
surface areas, number of municipalities and number of commuters. 

The last problem we solved is the one related to calibration. Applying the 
model with an extended job search base and an exponential deterrence func- 
tion, we found a constant equal to 1.94 • 10~^ is a perfect parameter value to 
generate commuting network of French administrative regions, whatever they 
are. However, we didn't test this result for other countries having different types 
of administrative regions. The robustness of this result to commuting network 
described at very different scales than the municipality one remain a question 
we want to address in the future. 
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Table 4. Description of the regions 



ID 


Region 


Number of 
municip. 

(region) 


Number of 
municip. 

(outside) 


Region area 

(W) 


Average 
municip. 

area (km ) 


Number of 
commuters 


FRl 


Auvergne 


1310 


3463 


26013 


19.86 


295776 


FR2 


Bretagne 


1269 


1447 


27208 


21.44 


653710 


FR3 


Ain 


419 


2809 


5762 


13.75 


162370 


FR4 


Alsace 


903 


3081 


8280 


9.17 


440961 


FR5 


Aquitainc 


2296 


2835 


41309 


17.99 


700452 


FR6 


Mayenne 


261 


3124 


5175 


19.83 


69915 


FR7 


Lozere 


185 


1859 


5167 


27.93 


12273 


FR8 


Poitou-Charente 


1464 


2467 


25810 


17.63 


375363 


FR9 


Centre 


1842 


4718 


39151 


21.25 


624693 


FRIO 


Midi-Pyrenee 


3020 


3845 


45348 


15.02 


o4o1d2 


FRll 


Limousin 


747 


3169 


16942 


22.68 


139481 


FR12 


Franche-Comte 


1786 


3317 


16202 


9.07 


268399 


FR13 


Hautc-Normandie 


1420 


3536 


12317 


8.67 


469335 


FR14 


Haute-Marne 


433 


3914 


6211 


14.34 


42690 


FR15 


Vosges 


515 


3808 


5874 


11.41 


92053 


FR16 


Lorraine 


2339 


3067 


23547 


10.07 


547457 


FRl 7 


Creuse 


260 


1814 


5565 


21.40 


23949 


FR18 


Langucdoc- 
Roussillon 


1545 


3046 


27367 


17.71 


409116 


FR19 


Charente-Maritime 


1948 


1983 


25606 


13.14 


375363 


FR20 


Haut-de-Seine 


36 


1245 


176 


4.89 


973173 


FR21 


Yveline 


262 


1543 


2284 


8.72 


618741 


FR22 


Val d'Oisc 


185 


1707 


1246 


6.74 


526600 


FR23 


Val de Marne 


47 


1234 


245 


5.21 


642092 


FR24 


Haut-Rhin 


377 


2283 


3525 


9.35 


183504 


FR25 


Tarn et Garonne 


195 


2338 


3718 


19.07 


41600 


FR26 


Pyrenee-Atlantique 


547 


449 


4116 


7.52 


65469 


FR27 


Alpes-Maritimes 


163 


353 


4299 


26.37 


163445 


FR28 


Loire 


327 


2788 


4781 


14.62 


178828 


FR29 


Territoire de Belfort 


102 


2031 


609 


5.97 


45185 


FR30 


Scinc-Saiiit-Denis 


40 


783 


236 


5.90 


655200 


FR31 


Essonne 


196 


1597 


1804 


9.20 


518321 


FR32 


Ardennes 


463 


2588 


5229 


11.29 


59963 


FR33 


Aube 


433 


2728 


6004 


13.87 


75561 


FR34 


Correze 


286 


2088 


5857 


20.48 


49815 



